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ABSTRACT 

Aims. We perform an advanced study of the fast magnetoacoustic sausage oscillations of coronal loops in the context 
of MHD coronal seismology to establish the dependence of the sausage mode period and cut-off wavenumber on the 
plasma-/} of the loop-filling plasma. A parametric study of the ratios for different harmonics of the mode is also carried 
out. 

Methods. Full magnetohydrodynamic numerical simulations were performed using Lare2d, simulating hot, dense loops 
in a magnetic slab environment. The symmetric Epstein profile and a simple step- function profile were both used to 
model the density structure of the simulated loops. Analytical expressions for the cut-off wavenumber and the harmonic 
ratio between the second longitudinal harmonic and the fundamental were also examined. 

Results. It was established that the period of the global sausage mode is only very weakly dependent on the value 
of the plasma-/? inside a coronal loop, which justifies the application of this model to hot flaring loops. The cut-off 
wavenumber k c for the global mode was found to be dependent on both internal and external values of the plasma-/?, 
again only weakly. By far the most important factor in this case was the value of the density contrast ratio between 
the loop and the surroundings. Finally, the deviation of the harmonic ratio P\j2Pi from the ideal non-dispersive case 
was shown to be considerable at low k, again strongly dependent on plasma density. Quantifying the behaviour of the 
cut-off wavenumber and the harmonic ratio has significant applications to the field of coronal seismology. 
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1. Introduction 

In recent years, improvements in instrumentation, data 
analysis techniques, and theory have led to confirmed ob- 
servations of magnetoacoustic waves in the solar corona, 
particularly the longitudinal, k ink, and sausage modes (e.g. 
iNakariakov fc Verwichtell2005l) . These waves attract atten- 
tion mainly because of the possible role they play in coronal 
heating and as natural probes of coronal plasma. The lat- 
ter concept gave rise to MHD coronal seismology, the novel 
approach to coronal plasma diagnostics by means of MHD 
waves. 

The sausage mode (also known as peristaltic or m = 
mode) is a symmetric perturbation of the loop minor ra- 
dius, in the form of periodic broadening and narrowing of 
a plasma tube (see Figure [T]). This mode, with typical pe- 
riods in the range 5-30 s, is believed to be detected in the 
microwave and har d X-ray emission associated with flar- 
ing solar loops (e.g. | Na kariakov ct al. 2003; Melnikov et al.l 
120051 : llnglis et al.ll2008 | ). and also in H„ emission from cool, 
post-flare loops (jSrivastava et al.|[2008F) . The sausage mode 
is essentially compressible, with the density perturbations 
in phase with the perturbations of the magnetic field and 
in anti-phase with the perturbations of the loop minor 
radius. This mode is a robust collective disturbance of a 
plasma structure, which is pract ically insensitive to small 
scale irregularities of the plasma (jPascoe et al.ll2007al) . The 
sausage mode can be manifest in both standing and trav- 
elling form. It is highly dispersive, with a pronounced de- 
pendence of the phase and group speeds upon the wave 
number. In the standing regime, this makes the spectrum 
of the resonant frequencies of different spatial harmon- 



ics non-equidistant. Long wavelength (in comparison with 
the minor radius of the oscillating loop) sausage modes 
are subject to a cut-off, in contrast with all other (kink 
and various fluting or ballooning) magnetoacoustic modes. 
Modes with wavelengths longer than the cut-off are leaky, 
and their phase s peed is slightly high er than the exter- 
nal Alfven speed (|Pascoe et al.|[2007b[ ). For shorter wave- 
lengths, sausage modes are trapped in the guiding struc- 
ture, and have phase speeds in the range between the in- 
ternal and external Alfven speeds. This cut-off wavelength 
is dependent on the density contrast between the loop and 
its surroundings. As flaring loops are usually dense, thick 
and short, both trapped and leaky regimes can occur. The 
main properties of the sausage mode are derived from the 
straight cylinder or slab model. Clearly, the slab and the 
cylindrical case s are n ot identical, as has been discu s sed in 
iKopvlova etaU (|2007l ). In particular. iTerradas et all (|2006t 
demonstrated the important difference in the wave leaking 
between these geometries. However, the difference in geom- 
etry should not be the principal difference in the resonant 
periods of sausage modes, as the behaviour of this mode is 
known to be almost identical in the slab and cylinder ge- 
ometries in the trapped regime. Future study of the sausage 
mode of a plasma cylinder will provide a final answer, but 
is out of the scope of this paper. 

Because of its compressibility, the sausage mode can be 
readily seen in microwave, EUV and soft X-ray bands, pro- 
vided sufficient time resolution is available. Also, this mode 
is capable of modulating the population of nonthermal elec- 
trons via changes in the magnetic mirror ratio in the flar- 
ing loop, leading to periodic precipitation of electrons to 
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the loop footpoint. This causes modulation of the emission 
in hard X-ray and white light because of bremsstrahlung 
(|Zaitsev fc Stepanovlll982D . 

The seismological potential of the sausage mode is con- 
nected with estimating the magnetic field outs ide the os- 
cillating plasma structure (jPascoe et al.ll2007bt ). and with 
determining the transverse profile of the plasma density in 
the flaring loop (jNakariakov et all [2004). Also, in the case 
when the microwave spectrum is generated by the gyrosyn- 
chrotron mechanism and is spatially resolved, this mode 
provides a promising basis for combined MHD-microwave 
diagnostics of flaring plasmas. As the sausage mode often 
causes pronounced modulation of the emission light curve, 
it also o pens up possibilit ies for the diagnostics of stellar 
coronae (|Nakariakovll2008D . 

One popular topic of MHD coronal seismology is the use 
of the ratios of global (or fundamental) modes and their sec- 
ond harmonics (variously referred to as the Pi/2P2-ratio or 
the Pi /Pi ratio) for determining the longitu dinal profile of 
the pl asma parameters in the oscillating loop. lAndries et ail 
(|2005f ) noticed that for the kink mode the discrepancy be- 
tween the observed P1/2P2 ratio and that theoretically cal- 
culated with the use of the straight cylinder model could 
be attributed to the variation of the plasma density along 
the loop, e.g. because of gravitational stratification. There 
have been several attempts to account for this stratification 
in the theoretical model (see iDvmova fc Rudermanl 12 00 



Donnelly eraHl2006l: iMcEwan et al.ll2008t IVerth fc Erdelvil 
2008) , and also to extend t his study to longitudinal modes 
(e.g. lDiaz fe Rober ts 2006). One of the difficulties with the 
full scale implementation of the Pi/2P2-based technique in 
MHD coronal seismology is that the higher harmonics of 
the kink mode are rarely detected because of insufficient 
temporal resolution of solar EUV imagers. 

As the sausage mode is usually observed with the instru- 
ments which have very high time resolution (e.g. 0.1 or 1 s 
in the case of the Nobeyama Radioheliogr aph) , it is possible 
to de t ect several harmonics of this mode (jNakariakov et all 
I2003t iMelnikov et al.l [20051) . However, as in contrast with 
the kink mode the sausage mode is highly dispersive, the 
identification of the sausage mode harmonic number is a 
non-trivial task. This is partly due to a significant depar- 
ture of the Pi/2P2-ratio from unity because of the disper- 
sive modification of the phase speed even in the absence 
of longitudinal stratification. It becomes especially difficult 
when spatial resolution is not available. Hence there is need 
for detailed theoretical modelling of various regimes of the 
sausage mode dynamics. 

It is not likely that in flaring loops the sausage mode 
is affected by gravitational stratification, as in such loops 
the density scale height is large because of the high tem- 
perature. The effect of the variable cross-section of the flar- 
ing loop on the resonant periods of the standing sau sage 
mode was recently considered by iPascoe et al.l ([2009D . It 
was established that the sausage mode continued to be 
supported despite variations in the loop cross-section. It 
was also shown that the characteristic periods of the global 
mode and higher harmonics were effected, but weakly so. 
A similarly small variation in parameters was found when 
inves tigating the effects of twist in loops ([Erdelvi fe Fedunl 
l2006h . However, the high temperatures and densities of flar- 
ing lo ops also lead to t he increase in the plasma-/? inside the 
loops (|Shibasakill200l . This can cause a significant devia- 
tion of the resonant periods and cut-offs from the results ob- 
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(a) The global sausage mode (b) The second (longitudinal) 

harmonic 



Fig. 1. Illustrations of a) The global sausage mode and b) 
The second (longitudinal) harmonic of the sausage mode, 
shown as filled density contours. In both cases the pertur- 
bation is symmetric about the loop axis, which remains 
unperturbed. 



tained i n the zero-/3 approxim ation studied in lCooper et al.l 
(1200.1 ): IPascoe et al.l (t2P^07bl) . 

The purpose of this study is to establish the effect of 
finite plasma- /3 on the resonant periods of the sausage mode 
by full numerical MHD modelling, and also to investigate 
the behaviour of the harmonic ratio for the sausage mode. 
The paper is organised as follows: in Section 2 we describe 
the standard theoretical approach to modelling wave modes 
in a magnetic slab geometry. Section 3 details the numerical 
methods used for the modelling in this paper. In Section 4 
we present our results, which are then discussed in more 
detail in Section 5. 



2. Sausage modes in a magnetic slab 

Following the standard approach, we model magnetoacous- 
tic oscillations in a coronal loop as perturbations of a field 
aligned enhancement in plasma density. The magnetic field 
is taken to be in the y-direction. The density profile is uni- 
form in the y-direction. In t he transverse direction the den - 
sity is given by the profile (jNakariakov fc R oberts 1995b) 



p(x) = (p - p e ) sech 2 



(1) 



where p e is the external density far from the loop, po is the 
density at the loop apex and a is a characteristic length 
governing the width of the loop. Here, n = 1 corresponds 
to the symmetric Epstein profile (see Figure [2]) which al- 
lows for analytical treatment. Increasing values of n lead 
to a steeper profile, meaning that with a large enough ex- 
ponent (n —¥ 00) the profile given in Equation (JlJ can be 
considered as an approximate step function, also with a 
known analytical solutio n for fast magnetoacoustic modes 
([Edwin fc Robertslll982D . Most importantly, the use of this 
profile avoids steep gradients that may lead to strong arti- 
ficial shocks in numerical simulations. 

The equilibrium condition is the total pressure balance 
between th e internal and external plasmas. This may be 
written as (|Edwin fc Robertslll982D 



p 1= 2C s 2 0+7 Ci 
Po 



(2) 
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where 7 = | is the adiabatic index, Cao and Ca c are the 
internal and external Alfven speeds, and C s o and C se are 
the internal and external sound speeds. 

For a symmetric Epstein profile, the solutions describing 
the perturbation of the transverse velocity component in 
the sausage mode are given by 



U(x) 



sinhfz/a) 
cosh A (a;/a) 



where 



(cic-c P 2 ) 1/2 fc 

Cab 



A = 



1, 



(3) 



(4) 



where k is the wavenumber and t he phase speed CV> = u/k 
is in the range C A0 < C p < C Ac ([Cooper et alJl2003h . 

Since C p is typically a large fraction of Cxe (for trapped 
modes) A is usually of order 1. The phase speed C p can be 
worked out analytically from the dispersion relation 



r 2 

°A0 



2 

ka 



3 

Cac 



r<2 en 



(5) 



This equation allows for the exact analytical solution for 
the phase speed C p . Unfortunately, this solution is only 
valid in the case of zero-/3 plasmas. 

For a realistic flaring loop the situation is likely to be 
more complex. A perfect pressure balance, for example, 
may not hold. Nor are the background plasma parameters 
likely to be uniform, or the amplitudes of oscillation suffi- 
ciently small to be considered linear. However, the global 
sausage mode is a very robust perturbation, and it remains 
desirable to examine its properties in a more idealised case 
before considering higher order effects. It is also likely that 
the choice of geometry is relatively unimportant for the 
MHD sausage mode. Part of the reason for this is its sim- 
plicity; a symmetric density perturbation about the loop 
axis independent of the azimuthal angle. This is in con- 
trast to other MHD modes. In the following analysis, we 
simulate sausage modes using both a step-function density 
profile and a classical Epstein profile, noting the rationale 
in each case. 



3. Numerical model 

The numerical simulations of MHD sausage modes pre- 
sented here were p erformed in slab geometry using Lare2d 
(jArber et al.ll2001f) . The code solves the MHD equations in 
normalised form. The defined normalisation constants used 
in this paper are: 

L Q = 1 Mm, 
B = 1 x 1CT 3 T(10 G), 
po = 7.95 x 1(T 13 kgixT 3 . 
This leads to the additional quantities: 

t = 1 s, 

v = 1 x 10 6 ms~\ 
T = 6 x 10 7 K. 
The normalised MHD equations may then be written as 



D P 

Dt H 



Dv 



p— = {VxB)xB- VP, 



DB 

~Dt 



(B ■ V) v — B (V ■ v) 



De 



(0) 



(7) 



(8) 



(9) 



where v is the plasma velocity, P is the gas pressure, B is 
the magnetic field, p is the mass density and e = Pj p(j — 1) 
is the internal energy density. 

An initial structure in equilibrium was set up as de- 
scribed in Section 2. In order to control the value of /?, 
the magnetic field was reduced inside the loop as shown in 
Figure^] All of the values in Figure [5] are normalised. The 
internal density of the loop, po, is also fixed. The internal 
gas pressure is left as a free parameter in order to maintain 
total pressure balance. As a result the plasma temperature 
varies with the internal gas pressure. Consequently, the sim- 
ulated plasma structures are both hotter and denser than 
the surroundings, consistent with the known properties of 
physical coronal loops. 

A convenient way to modify the equilibrium is to ad- 
just the profile of the magnetic field, keeping all other pa- 
rameters the same. This is the approach that will be used 
throughout this paper. This would allow us to consider the 
important case of flaring loops filled in with a hot and 
dense, high-/? plasma. 

In the following, the values of plasma-/? in the centre 
of the loop and at infinity are referred to as /3o and /3 e , 
respectively. To simulate a sausage mode we set up an initial 
transverse velocity perturbation of the form 



v x = AU(x)sin(ky), 



(10) 



where A is an amplitude and x and y are the coordinates 
across and along the slab respectively. This velocity dis- 
tribution disturbs the loop at its centre and ensures that 
a global standing mode is developed. To avoid nonlinear 
effects the amplitude of the perturbation must be small 
compared to the equilibrium parameters. We choose A = 
0.002 in all cases, resulting in perturbations which are very 
small fractions of the equilibrium magnetic field and den- 
sity (typically ~ 0.1% of the internal density for example). 
As described in Section 2, the magnetic field is aligned in 
the y-direction, as is the loop itself. Reflective boundaries 
were applied in the y-direction, while open boundaries were 
applied in the transverse direction. 

Theoretically, any excitation of a mode in this 
way would lead to the excitation of higher harmonics 
(jTerradas et al.ll2005D . However, the excitation we choose 
is always sufficiently close to the eigenmode that the am- 
plitude of these harmonics is very small with respect to 
the fundamental oscillation, and they are quickly damped. 
Thus, by running the simulation for suffici ent time, the 
domin ance of the desired mode is ensured (jPascoe et al.l 
l2007bT) . 

To measure the periods of oscillation, the perturbation 
of the plasma parameters at the relevant anti-node are stud- 
ied. For the global sausage mode this is at the loop centre. 
The period is determined by applying a Fourier transform 
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Fig. 2. Top panel: plot of the normalised initial parameters 
as a function of the transverse coordinate. The solid line is 
the density, the dashed line is the magnetic field B y , and 
the dot-dashed line is the temperature. The temperature 
has been multiplied by a factor of 5 for clarity. Bottom 
panel: the j3 profile resulting from the parameters shown in 
the top panel. 



Fig. 3. Period of the global sausage mode as a function of 
loop length L, for a number of f3 values. The stars corre- 
spond to /3 =0.005, the diamonds to /?o=0.08 and the tri- 
angles to /?o=0.82. The solid line corresponds to 2L/Cas, 
while the dashed line is the analytical solution to the dis- 
persion relation for an Epstein density profile. The density 
contrast ratio is 10. 
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to the time series data from this point, thus the uncertainty 
in the measurement is connected with the half-width of the 
peak in the Fourier spectrum. 



Fig. 4. The sausage mode period dependence on /3 for a 
trapped mode in a loop of length L=15 Mm. The density 
contrast ratio is 10. 



4. Results 

4.1. The period of the sausage mode 

According to Nakaria kov et all (|2003|) . the period of the 
global sausage mode near the wavenumber cut-off in the 
trapped regime is given by 



P 



2L 



(11) 



where L is the loop length and C p is the phase speed. More 
recently, this relation w as shown to be a val id approxima- 
tion in the leaky regime (iPascoe et a l. 2007b) and in a loop 



with variable cross-section Pascoe et al. (|2009l ). However, 



in both cases this relation was tested for very low values 
of (3q. It is unclear whether or not a significant value of /?o 
would have a strong effect on the period of th e mode. As a 
first step therefore, we reproduce the result of Pascoe et al. 
(|2007bl) showing the dependence of the period on the loop 
length L, but this time for a variety of (3q values. This is 



achieved by simulating sausage oscillations in a loop with 
an Epstein density profile and a density contrast ratio of 
10. 

The results shown in Figure El are in excelle nt agreement 
with those presented in IPascoe et al.l (|2007bT i . Above ap- 
proximately 20 Mm the mode becomes leaky and the phase 
speed is slightly higher than the external Alfven speed. 
Furthermore, the near-independence of the oscillation pe- 
riod on Pq is clear. As a further illustration, we investi- 
gate the period dependence on fio at a fixed length L. The 
length of the loop was set at L=15 Mm - corresponding to 
trapped modes - and the value of /?o was varied as described 
in Section 3. /3 e was kept constant and small throughout. 
The result is displayed in Figure 01 

The period of the sausage mode should also depend on 
the loop width a. This is due to dispersion, since a change 
in a essentially corresponds to movement along the disper- 
sion curve derived from Equation [5] and hence a change in 
phase speed C p . To illustrate this, we performed numeri- 
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Fig. 5. The sausage mode period dependence on width a for 
a loop of fixed length L=15 Mm, where the density contrast 
ratio is 10. The solid line corresponds to the solution of the 
dispersion relation for an Epstein profile. 



cal simulations of a sausage mode in a loop of fixed length 
L=15 Mm while varying the width of the Epstein profile. 
The results are illustrated in Figure[5l As expected, the pe- 
riod increases as a function of a within the trapped regime, 
whereas for leaky modes - corresponding to a < 2.2 Mm - 
the period tends to a constant value. This is a consequence 
of the phase speed remaining close to the external Alfven 
speed in this regime. 

It is clear from these results that, although there is a 
slight increase in the period with respect to /3rj, the variation 
is very small (less than 5%) compared to the effect of other 
parameters such as loop length. Thus we may say that the 
period of the sausage mode is very stable with respect to 
the value of /?o- 

4.2. The second longitudinal harmonic 

An often-neglected aspect of sausage mode behaviour is 
the relationship between the periods of separate harmon- 
ics. For dispersionless standing modes in a resonator, the 
frequency of the AT-th harmonic is simply N times that 
of the fundamental. Being strongly dispersive, the sausage 
mode harmonics deviate significantly from this pattern. 

The true behaviour of the harmonic ratio is explored 
for the fundamental and second harmonic modes. For a 
given value of k\ (the wavenumber corresponding to the 
fundamental mode), the period ratio P1/2P2 is measured 
from numerical simulations. By varying the loop length L, 
we control the wavenumber k\ and subsequently measure 
the period ratio as a function of k\ . Figure [6] displays the 
results. 

For non-dispersive waves, the expected ratio of P\ /2P2 
is 1. It is clear that the sausage mode harmonics do not 
follow this pattern. For large values of k\ corresponding to 
the uniform plasma (and hence dispersionless) limit, the 
ratio does indeed tend towards the ideal value. However, 
for smaller ki, the period of the second harmonic be- 
comes a much greater fraction of the fundamental period, 
as illustrated in both panels of Figure El The top panel 
shows the quantitative difference between simulations of 
an Epstein profile (diamonds) and a step-function pro- 
file (stars). Although this deviation is noticeable at low 
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Fig. 6. Top panel: The ratio P\j2Pi observed via full nu- 
merical simulation for the sausage mode as a function of 
the longitudinal wave number k\ for both the step-function 
(stars) and Epstein profile (diamonds) geometries. The den- 
sity contrast ratio is 20 and /3o=0.08. The solid line is the 
theoretical curve derived from the dispersion relation for 
the Epstein profile, also using a density contrast ratio of 20. 
Bottom panel: Illustration of the dependence of the analyt- 
ical solution on the density contrast ratio. The solution is 
shown for contrast ratios of 5 (upper solid line), 10 (dashed 
line), 20 (solid line), 50 (dot-dashed line), 100 (dotted line) 
and 1000 (triple-dot-dashed line). 



k\ values, the overall behaviour remains the same, as ex- 
pected. The solid line shows the theoretical curve for an 
Epstein profile derived from the dispersion relation given in 
Equation (|5|), corresponding to the zero-/3 case. The agree- 
ment between the theoretical and numerical results is im- 
pressive. The bottom panel of Figure [6] shows the same the- 
oretical curve for a range of density contrast ratios. Clearly 
there is a strong effect on the behaviour of the period ra- 
tio; the denser the loop, the greater the deviation from the 
ideal, non-dispersive case. However, the range for the pe- 
riod ratio in all cases lies between 0.5 and 1. The turnover 
of the period ratio curve seen at low k\ is a product of being 
very close to the cut-off wavenumber k c . 

In coronal seismology, the ratios of observed oscillatory 
periods - say, in a coronal loop or solar flare - are often used 
to make judgements about the underlying MHD modes. 
Figure [6] provides a powerful illustration of the dispersive 
nature of sausage mode harmonics, and that for entirely 
realistic values of k it is not expected that the harmonic 
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ratio be close to 1 regardless of any non-uniformity in the 
longitudinal direction. 

4.3. The cut-off wavenumber 



According to iNakariakov fc Roberts! (|l995al ). the cutoff 
wavenumber for a sausage mode in a magnetic slab with 
a step function profile of finite-/3 plasma is given by 



k c a = — 
2 



\ 



. Ci +Cjf 



c- 

A 



i 



s)( 



r- 



An analytical expression for the wavenumber cut-off is 
not readily obtainable in the finite-/? case for the symmet- 
ric Epstein profile. The available expression (Pascoe et al. 
2007) is only valid when /3=0. Therefore for consistency all 
analysis of the cut-off wavenumber is carried out in step- 
function geometry. 

In order to investigate the dependence of this cut-off 
wavelength on the plasma- /3, we write the sound speed as 
C S 2 /C A = 7/3/2. Equation [12] then becomes 



k r a = — 



\ 



<■- 



2 



r 2 

A 



1 - 



\ 2 cTo) 



(13) 



(14) 



Using the equation for magnetostatic pressure balance 
(Equation [5]), we can eliminate Cas/Cao and write the ex- 
pression for the cut-off wavelength in terms of the density 
contrast £ = po/p c and the plasma- (3. We find 



r* 2 

°Ac 
°A0 



C(l + /3o)~/3c 
Inserting this expression into Equation 1141 we obtain 

(CA) - & + Cf 




- fl, + C - 1) ((C - i)/3 - /?c + C) ' 

This dependence is shown in Figure [7J For any value of 
/3 e , k c a is seen to decrease as a function of increasing /3q. 
Realistically, in a coronal plasma the external j3 e is likely 
to be small. Thus we will concentrate on this region of the 
surface. 

By comparison, we numerically evaluate Eauation ll2l for 
increasing values of j3o, while maintaining a fixed and low 
j3 e . We maintain fixed external values of p and B. Inside the 
loop, the value of B is decreased, increasing the value of /3q. 
For the structure to remain stable, the internal gas pressure 
must be left a free parameter to maintain total pressure 
balance. As before, the temperature is allowed to vary with 
the gas pressure, ensuring that the simulated loops are hot 
and dense. The results are shown in Figure [51 where the 
solid line corresponds to a density contrast ratio of 10. 

In both cases, the value of the cut-off wavenumber 
clearly decreases while /3q is increased from zero to one. 




Fig. 7. Surface showing the dependence of the cut-off 
wavenumber on /3q and /3 e . The density contrast ratio is 
10. 




0.0 



1 .0 



Fig. 8. Plot showing the dependence of k c a on 0q for a fixed 
value of /3 e , where /3 e = 0.01. The density contrast ratios 
used are 10 (solid line), 15 (dotted line), 25 (dashed line) 
and 50 (dot-dashed line) . The overlayed datapoints are the 
results from full numerical simulations. 



Figure [8] shows the same result for a number of different 
density contrast ratios, illustrating that the density has a 
strong impact on the confinement of the sausage mode. 
Again though for any starting - and fixed - value of the 
density contrast, the cut-off wavenumber will decrease as a 
function of /3q. 

To further corroborate these results, we perform full nu- 
merical simulations of the global sausage mode for varying 
values of /3q. This is achieved by controlling the internal 
magnetic field as explained in Section 2. 

To measure the cut-off wavenumber itself, we simulate 
the global mode for varying L (and hence k) and exam- 
ine the time series of the loop anti-node as described in 
Section 3. Here the criterion for a mode to be leaky is that 
the amplitude of the oscillation inside the loop is damped. 
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When far into the leaky regi me, the damping effe ct is strong 
and evident (see for example iPascoe et al.l2007bT ). However, 
closer to the trapped regime it can be difficult to ascertain 
whether significant damping is taking place. This effect 
is mitigated by running the simulations for a sufficiently 
long time. However an exact determination of the cut-off 
wavenumber k c remains difficult. 

The results of these MHD simulations are overlayed in 
Figure [5] The simulations show the same general trend as 
the analytical and numerical solutions of Equation 1121 al- 
though the agreement is not exact. The curve of the sim- 
ulated result is noticeably shallower than the theoretical 
predictions. In part this is likely due to the difficulty in de- 
termining exactly the transition between the trapped and 
leaky regimes, a problem that is reflected in the size of the 
error bars shown for the numerical results. 



5. Conclusions 

Using numerical simulation and analytical study, we have 
carried out detailed analysis encompassing various aspects 
of MHD sausage mode behaviour. The key findings may be 
summarised as follows: 



1. 



In the finite-/? regime, the period of the global sausage 
mode was found to be determined by the length 
of the loop L, which c onfirmed previous findin gs of 
Nakariakov et all (l2003h . lAschwanden et al.l (|2004f ) and 



Pascoe et all (|2007bh . 



2. In the trapped regime the period is also dependent on 
the loop width a due to dispersion, but becomes weakly 
dependent in the leaky regime. 

3. The period of the global sausage mode is not signifi- 
cantly affected by finite values of /3q, with the variation 
remaining less than 5% for < (3 < 1. 

4. The ratio between the second longitudinal harmonic and 
the fundamental, P1/2P2, varies between approximately 
0.5 and 1 for 'reasonable' values of k\ and p, and is 
strongly dependent on both k\ and p. 

5. The cut-off wavenumber for the global mode is a func- 
tion of (3q, (3 e and the density contrast ratio £. For a 
given density and /3 e , the k c a is a decreasing function 
with respect to /?o- 

Analysis of the cut-off wavenumber k c reveals the lim- 
iting factors dictating whether an MHD mode remains 
trapped or becomes leaky. The unambiguous result is that 
the density contrast between the loop itself and the sur- 
roundings is the most important factor. For extremely 
dense loops, a sausage mode may remain trapped even for 
very long loop lengths. Conversely, for loops of modest den- 
sity, it is clear that trapped modes may only be supported 
for comparatively short lengths. The values of f3 e and /?o, 
although also factors, are secondary ones. It was shown in 
Section 14.31 that k c decreases as a function of /?o when /3 e 
and C are treated as constants. In more physical terms, 
this means that the maximum length supporting trapped 
modes, L Cl is increased. However, this variation is clearly 
small when compared to the effect of the density contrast 
ratio Q (see Figure [5]). 

These results also have immediate practical applica- 
tions in the context of MHD coronal seismology. For multi- 
periodic oscillations, an understanding of the relationship 



between sausage mode harmonics enables informed judge- 
ments about the na ture of the observed mo de. 

For example, in ISrivastava et al.l |2008) a period ratio 
of Pi/2P 2 Kt 0.83 was reported, with the favoured inter- 
pretation being the fundamental and second harmonic of 
the sausage mode. From their estimates of the loop length 
(L w 100 Mm) and our independent estimate of the loop 
width (a«8 Mm) , it is possible to examine this conclusion. 
Given the values of L and a, we estimate that k\a w 0.25. 
It is clear from Figure H2 that, for the period ratio reported 
this is very much in the leaky regime for the sausage m ode. 
However, the data presented in lSrivastava et aTl (|2008l ) does 
not show significant evidence of damping. Therefore, al- 
though it must be pointed out that the uncertainty in k\a 
is considerable, it seems that the observed P1/2P2 ratio is 
an indication that another mode may be responsible. 

We apply the same analysis to the event analysed in 
iNakariakov et all (|2003|) . where a fundamental and possi- 
ble second harmonic was reported. In this case the observed 
period ratio is P1/2P2 = 0.82 ± 0.15. Given the published 
estimates of the loop length (L w 25 Mm) and width (a » 
6 Mm) we find that in this case k\a k, 0.75. Again, we 
can compare this with the results shown in Figure [51 In 
this case, the observed period ratio and wavenumber are 
broadly consistent with a trapped fundamental and sec- 
ond harmonic of the sausage mode, provided the density 
contrast ratio is not too large, say 5 or perhaps even 10. 
Although this is not particularly large for a hot flaring loop, 
it is certainly possible and consistent with a trapped mode 
regime. Another possibility to consider is that the observed 
oscillations are actually triggered by a nearby non-flaring 
loop, which might have a lower density c ontrast than the 
flaring loop itself (|Nakariakov et ai1l2006t ). These two ex- 
amples show how the results obtained here can be readily 
applied to observational examples of multi-period QPP. 

Nevertheless, significant scope for improvement re- 
mains, both in observational and theoretical terms. In gen- 
eral, the uncertainties in both P1/2P2 and k\a are sub- 
stantial in observational studies of QPP. Additionally, as 
has been shown in Figure H2 the geometry of the theoret- 
ical model has a bearing on the quantitative values of the 
period ratio. Potentially important effects such as loop cur- 
vature and gravity have also been neglected in this study. 
Understanding the impact of these effects may prove an 
important next step in quantifying the characteristics of 
sausage modes. 
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